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THE ANALYSIS EXPERIMENTS 
INVOLVING THE KINETICS PILES WITH 


Rumsey? 


ABSTRACT 


The theory pile kinetics used analyze measurements the mean life- 
time and migration length neutrons. The purpose show what effects 
need taken into account the analysis pile experiments rather than 
give precise experimental verification the theory. shown that the bare- 
pile one-group model not adequate. two-group model which takes into 
account the presence the reflector used obtain more accurate analysis. 
The orders magnitude are illustrated deriving formulas for specific pile— 
the Chalk River ZEEP. found that the presence the reflector results 
increase the mean lifetime about 20%. Increasing the multiplication 
constant given amount produces bigger value effective multiplication 
constant for the pile with reflector than without, but owing the larger 
value the pile with reflector diverges more slowly than the bare pile. The 
bare-pile model overestimates about 10% the change required 
re-establish stability after given change The migration length over- 
estimated about 25%. The more elaborate theory agrees well with various 
measured values, but this partly fortuitous because the structure the ZEEP 
rather complex and not adequately represented the theory. 


INTRODUCTION 


Various experiments the ZEEP have shown that the simple one-group 
model based the assumption bare pile (3, for the analysis pile 
kinetics not consistent with the actual behavior ZEEP. This not sur- 
prising because the one-group model applies strictly bare pile for which the 
reproduction constant very slightly greater than unity. The experiments 
Fergusson and Gilbert (2) bring out one aspect the matter. These investi- 
gators measured the effect the pile power absorber whose rate 
absorption was made oscillate sinusoidally with periods from 0.5 500 sec. 
They found that the mean lifetime neutron derived from their results was 
significantly greater than the value calculated from cross sections, etc. the 
basis the bare-pile model. The experiments Bayly (1) bring out another 
aspect. Here, the time variation pile power was measured after the depth 
heavy water, denoted had been changed known amount from the 
value which gave steady power. found that the migration length neutrons 
derived from the experiment significantly greater than the value obtained 


other means, again assuming that bare-pile model applicable. 
received January 25, 1954. 
Contribution from Atomic Energy Canada Limited, Chalk River, Ont. Issued A.E.C.L. 
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Our aim is, then, analyze experiments the following type. Suppose that 
the pile running steady power and that suddenly increased 
amount 6h. result the pile power will start rise. wish determine 
the relation between the way which the power changes with the time and 6h. 

order give point the method described here, will review the pro- 
cedure which has been adopted such cases. argued that increasing 
above the critical height equivalent increasing the multiplication 
constant defined with respect lattice extending infinity. The increase 
k,, calculated from formula the form 


some more refined form where represents the Laplacian, which deter- 
mined from and the migration length the fissile medium. assumed 
that the amount which k,, would have reduced order re-establish 
steady power equal the value where represents the effective 
multiplication constant height above critical height. This one assump- 
tion that introduces some error. Again assumed that the mean lifetime 
which goes with equal the mean lifetime neutron the infinite 
medium, that is, the average time between the birth and death neutron 
spent all its life the heavy water uranium lattice. This procedure 
accurate for infinite pile but only approximately valid for piles finite size. 
For instance, does not take into account the fact that the average lifetime 
neutron somewhat lengthened because significant number neutrons 
pass through the reflector before being finally captured the lattice, and, 
owing the small capture cross section the reflector, neutron can spend 
relatively long time over this process. 


METHOD ANALYSIS 


The difficulty with the application the simple theory appears mainly 
due the fact that quantities like and are not clearly defined for pile 
such the ZEEP, which neither nor infinite: consists tank 
containing the heavy water uranium lattice and graphite reflector. can 
shown (6), the basis two-group model, that the simple equations 
can still retained for such cases, provided define the effective multiplica- 
tion constant and mean lifetime follows. 

The mean lifetime used the average time between the birth and death 
neutron the tank. This average time will average over all possible 
ways which neutron starting out from fission somewhere the tank 
finally gets captured thermal neutron the tank, and obviously 
influenced the presence the reflector. give some idea the orders 
magnitude involved, note that the mean lifetime neutron which spends 
all its life the tank about 0.7 msec., whereas the mean lifetime 
neutron which spends all its life the reflector about msec., and the 
actual mean lifetime, which have defined above some average between 
these two values. 

Let denote the lifetime neutron which spends all its life one 
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medium the Greek letter Thus property the medium; the case 
monokinetic neutrons the capture mean free path divided the 
neutron velocity. 

The effective multiplication constant which has used conjunction 
with the mean lifetime defined times the probability that neutron 
that born the tank will die the tank. Since this probability found 
taking the average over all possible ways which the process takes place, 
also depends upon the presence the reflector. Note that neither nor 
dependent any way the fact that there may delayed neutron sources 
sources well the ordinary prompt fission sources. can make use 
this order simplify the derivation formula relating and the 
amount which would have changed order re-establish stability. 
consider pile which identical the actual pile all respects, except 
that all neutron sources are prompt sources; particular, the distribution 
neutron sources the same the actual and piles. follows that 
and are the same these two piles. 

order illustrate the method adopt very simple model. Suppose 
that the pile spherically that the diffusion coefficients are equal 
the tank and reflector, and that one-group approximation valid. 
specific, assume that Then the prompt pile diverges like exp where 


This follows immediately from the definitions and Now whether not 
the pile diverging, the neutron density solution equations like 


and the critical condition which always has satisfied 
[4] 


where the radius the tank. some means and are varied, but 
remains constant, then 


2 


now consider two cases—(1) the diverging prompt pile, and (2) the 
same pile which have reduced k,, say order re-establish 
stability. take for the values and [5] the difference between the 
values and the first and second cases. For the diverging pile 
readily verified that 


where equal the value the tank, and the migration length 
the tank, and 

where equal the value for the reflector and the migration 
length the reflector. For the static pile have 


| 
| 
| 
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and 
[9] 
From [8] and [6] obtain 


and from [7] and [9] obtain 


Note that the assumption equal diffusion coefficients the tank and 
reflector implies 


Substitution from [2], [10], and [11] [5] gives the following relation between 


which have denoted the steady state value order conform 
the usual notation. 

will noted that contrary the usual assumption, and k,, 
are not equal but differ amount which diminishes the size the pile 
increases. 

need hardly added that the change k,, required for the prompt pile 
identical that required for the actual pile since k,, does not depend 
delayed neutrons, etc. 

can now determine the value due raising above the critical 
value amount 6h, finding the amount which would have 
reduce order have stable pile height 6h. The actual evaluation 
this relation given below, where consider the specific case the ZEEP. 
The formula [13] not valid the ZEEP because the simplifying assumptions 
made order derive not hold. 

The method evaluating the mean lifetime can illustrated means 
the following simple model. Suppose that all neutrons the system have the 
same velocity. the total assembly neutrons the whole system (tank and 
reflector) let the number the tank and the number destined die 
the tank. (Note the difference—many neutrons the reflector may destined 
die the tank.) Then 


[14] 


where given equations [3] and the integral taken over the volume 
the tank. find the value i.e. the number neutrons destined die 
the tank, regard the actual neutron distribution distribution 
sources, and solve for the new distribution which would arise from this dis- 
tribution sources. other words, given 


[15] 
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taken over the volume the tank, where 


Equations [16] imply that the actual number neutrons given 
volume, then the number that are destined die that volume. 

Now can express terms these quantities noting that the rate 
which neutrons are dying the tank can represented either 
n/T. Equating these two expressions and substituting for and from [14] 
and [15] obtain the following formula for 


the integrals being taken over the volume the tank. 

The mean lifetime fundamental constant the system. can 
measured using the delayed neutron periods time standard. the 
experiments Gilbert and Fergusson, which have referred, measure- 
ment the phase difference between the and power oscillations function 

The quantity can measured analyzing the variation pile 
power with time. this means the value obtained relative 
where that fraction the total yield neutrons from fission that are 


delayed. 
THE VALIDITY THE ONE- AND TWO-GROUP REPRESENTATIONS 

The one- and two-group models are usually compared with reference 
static system. this basis the simple one-group model works quite well for, 
say, the ZEEP. However the approximations involved time dependent 
problem are not the same kind static problem. For instance the static 
solution depends the slowing down length, and thermal migration length, 
about equally, because most cases and are approximately equal. 
time dependent case the times taken for slowing down and thermal migration 
are important and since the slowing down time much smaller, the thermal 
migration relatively much more important the time dependent problem 
than the static problem. 

The different nature the approximations involved the static and kinetic 
problems brought out even the simple case bare pile. The values 
obtained for the basis one- and two-group models this case are sig- 
nificantly different. obtain: 


one-group 


probability escaping resonance capture uranium, 


where 
and 


the one-group the sum and thermal and and are the 
slowing down and thermal diffusion lengths. Putting and using the 
approximate relation 
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have 


k,, the ZEEP, there 10% difference, but what more impor- 
tant from our point view that the correction due the pile not being 
infinite (in which case would differs factor two between 
the two models. Thus clear that the one-group model inadequate for the 
purpose calculating the other hand the formula relating and 
the change has the same form for bare pile however many groups 
use, i.e. 


APPROXIMATIONS USED THE FORMULA FOR 


order reduce the algebra while retaining sufficiently accurate repre- 
sentation, adopt the following device for the calculation note that 
the time spent slowing down either the tank the reflector less, 
order magnitude, than the time during which neutron migrates with 
thermal velocities. Furthermore, the slowing down times tank and reflector 
are more nearly equal than the thermal lifetimes. can therefore put 


where the value corresponding slowing down the tank and the 
average time between birth thermal neutron the system (tank and 
reflector) and its death the tank. Thus calculated means the 
formula [17], i.e. 


where 
[18] x = L? V2X + Trt Pt in the tank, 


these equations the density thermal neutrons. 

determine therefore have find sufficiently accurate representa- 
tion the thermal neutron distribution the system. well known, the 
one-group model gives good representation the tank except points 
within migration length the tank wall; the reflector the form 
again agreement with measurement except points near the tank. However 
can shown that the one-group model underestimates the number ther- 
mal neutrons the reflector relative the number the tank, factor, 
say which about equal two, compared with two-group model 
models more than two groups. order evade the turgid algebra asso- 
ciated with the two-group calculation, therefore take being given 
the one-group density the tank and times the one-group density the 
reflector. This then gives sufficiently accurate representation the relative 
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numbers tank and reflector and also the spatial distribution tank and 
reflector except for small region near the tank wall. 

The factor depends the relaxation lengths tank and reflector. 
not much greater than unity practically independent the shape the 
boundary between tank and reflector, and can therefore evaluated solving 
for the simplest possible geometry, e.g. for plane slabs spheres. Assuming 
equal diffusion coefficients tank and reflector the two-group model gives 
value for which can simplified the following approximate form: 


which accurate within few per cent 1.2; the thermal diffusion 
length the tank. 

APPROXIMATIONS USED DERIVING FORMULA FOR 


calculate the formula relating and the equivalent change 
the two-group critical condition should used place equation [4] but 
this impracticable. Since the critical condition depends more over-all 
effects rather than individual neutron groups, the one-group model 
satisfactory for this part the calculation. Note that even modify the 
one-group model means the factor introduced above, the critical 
condition remains unchanged. This procedure does not imply that 
‘insensitive the way which the and thermal neutrons are distributed, 
but rather that one-group model good enough for the purpose establishing 
the critical condition, especially since shall fit the given pile geometry. 

The next step deriving the formula for substitute values for 
and the differential relation obtained from the critical condition, 
cf. equation [8]. equations [6] and [7] show, and depend time 
constants, e.g. wr, [6], and consequently must cautious about using 
the one-group model here, because the inaccuracies mentioned earlier. 
Actually not difficult evaluate and means the two-group 
model. Adopting our previous notation which unprimed quantities refer 
the diverging prompt pile and primed quantities the prompt pile which 
the two-group model gives 

and place equation [8], 


where 
are the fast and thermal diffusion lengths the tank, 
the fast and thermal values the tank, 


Equation [10] then replaced 
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The two-group model gives the same relation [11] for but with slightly 
different interpretation: 
[22 2A 6A = 


where the thermal diffusion length and the thermal the reflector. 

determine our procedure then use the one-group critical 
condition but interpret the parameters and therein according the 
two-group model. 

CALCULATION FOR THE ZEEP 

The ZEEP described Reference (4). The ZEEP tank cylindrical. 
assume that the reflector infinite and covers the curved sides only the 
cylinder. Actually there some reflector the bottom the ZEEP tank, 
but this taken into account some extent fitting the constants give 
the measured values neutron distribution. 

The details the calculation are given the Appendix. obtain the 
following formula for assuming equal thermal diffusion coefficients tank 
and reflector: 


where 

and are the fast and thermal values the tank, 
Laplacian, 

thermal diffusion length the tank, 
thermal diffusion length the reflector, 

the factor defined above equation [20], 

tank radius, 

effective tank height, 


factor p), 
and are Bessel functions imaginary argument ipa, and 


t 


and are Bessel functions imaginary argument 


must emphasized that this formula should used only for assessing 
orders magnitude; any agreement with experimental results will 
coineidence. This partly due the fact that have not represented the 
ZEEP properly—the heavy water reflector the bottom the tank has not 
been allowed for all well—and partly due inadequate information about 
the The constants are assumed have the following values: 


7.9 inverse meters”? (measured value—Johns and Sargent (5)), 


meter, 


1.5 meters, 
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thermal diffusion length reflector 0.5 meters (measured value), 


125 based the assumption that 1.2 (by 
lattice calculation) and that 


115 cm.? The slowing down length, the lattice tied 
measured value the indium resonance 110 cm.? for pure heavy water and 
cannot far out. then fixed through the measured Laplacian. 


The controlling factors and cannot estimated all accurately. 
have taken 


(fast) 


where the transport mean free path and the velocity thermal neutrons 


(3). 


(fast) 
gives 0.05 msec. 
2 


(thermal) 2.4 cm. 
giving 
0.71 msec. 


Our formula then becomes 


The first two terms are what the value would there were reflector. 
see that the reflector constitutes correction the order 20%. 
use the above values for and then 


0.84 msec., 


which happens agree fairly well with the value 0.86 msec. obtained 
Fergusson and Gilbert (2). 


ASSOCIATED CALCULATION 
The assumptions listed above give the following formula for (see 
Appendix): 
Ree k 1 9 9 


which the and Bessel functions have arguments and respectively. 
Substitution the numerical values gives 


The expression the left what the time constant the prompt pile would 
the pile were infinite, the expression the right outside the brackets 
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the actual time constant the prompt pile. The second term brackets 
due the finite pile size and the third due the reflector. Substitution 
the calculated values and gives 
interesting compare this with what the result would the ZEEP had 
reflector. this case 


Increasing given amount thus produces bigger value for 
the pile with reflector than without, but owing the much larger value 
the former case, the pile with reflector diverges more slowly than the bare 
pile. worth noting that are given then the bare-pile model 
overestimates k,, about 10%. 

The relation between k,, and 6h, the increase height over the critical 
value, needs careful examination for the ZEEP. For this purpose seems 
desirable make specific allowance for the heavy water reflector the bottom 
the tank. consider the theoretical aspect order bring out the rele- 
vant details but here again theory cannot make very accurate prediction. 
does however suggest way which the desired relation can obtained 
experimentally. 

have find expression for where and are the Laplacians 
the same pile, which the unprimed case the actual ZEEP and the 
primed case what the ZEEP would reduced its k,, order 
make stable the increased value 

the cylindrical case, the radial and axial curvatures the density are 
specified and (see Appendix) where m?. assume that 
the radial conditions are the same any value (There insignificant 
error this assumption which due the fact that the height the reflector 
not kept line with the height core.) follows that 

illustrate the point suppose that the reflector the bottom the tank 
uniform thickness and characterized thermal diffusion length 
Let distance measured parallel the pile axis from the bottom the 
uranium rods, i.e. the reflector the bottom the tank extends from 

Adopting one-group model, the density has the form 

sin m(h the lattice, 


Since the diffusion coefficients are equal either side, the boundary con- 
dition has the form 


the height changed see that 
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where the ratio dp/dz the plane through the bottom the 
uranium rods and quantity that can determined experimentally. 


9 = 
The effective height given 
[27] 


and found analysis the vertical neutron distribution. put 
obtain the condition when there reflector below the rods. this case 
Hand 


Taking the following values 
150 cm.—effective height, 
120 cm.—height above bottom rods, 


obtain from [25] and [27] 
35cm. 


and equation [26] has the form 


compared with the relation [28] obtained ignoring the reflector: 


Ignoring the reflector therefore predicts value for due toa 
given that about too small. 
take given the two-group model, i.e. 


then 
Note that this differs from the one-group assumption 


about 10%; the one-group model predicts, from given and k,, 
value that too large about 10%. 
Summarizing the steps involved finding the relation between and 


have: 
(a) Connection between and k,, 
(c) Connection between and 6h, 


and each these steps there difference from the conventional approxi- 
mations the order 10%. 


4 
| | 
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terms our equations the steps are: 


[24] (a) ke 1 6k. Te E + K Ly + Vie (Ji/Jo)° + 1 ’ 


from which obtain the following formula for 


2 Qr2r7r 2 k. —1 Te 27 2 M’ (Ky, zy — 


illustrate the orders magnitude suppose that and are known, 
then the results obtained from the approximations listed below are: 
(i) Bare Pile—one-group model 
step overestimated about 10% 
step overestimated about 10% 
Total overestimated about 25% 


(ii) Bare Pile—two-group model 
step overestimated about 10% 
step significant error 
step (c) overestimated about 
Total overestimated about 15% 


(iii) Bare Pile—age velocity model 
This differs from (ii) one two per cent. 


The tendency for the estimate come down better models are 
used, cf. drop from 25% 15% between (i) and (ii). The age velocity model 
applied bare pile overestimates about 10%. 

Substitution the numerical values which have assumed gives 


De 


1.03 


had ignored the presence the reflectors would have obtained the 
erroneous result 


Formula [30] derived from two-group model. slight improvement 
accuracy can obtained replacing the expression for 
obtained from the age velocity model. 


on a7 
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APPENDIX 
The solution equations [3] written terms Bessel functions 


distance from axis, 
distance along axis, 
which satisfy the boundary conditions +3H,r 
Assuming equal diffusion coefficients tank and reflector and equating 
densities and currents the tank wall have 


[104] BKo(pa), 

[105] 

giving the critical condition 


Formula for 

explained the text the density thermal neutrons represented 
(equation 100) the tank and (equation 102) the reflector. Evi- 
dently the dependence the same for all densities that can 
omitted what follows. The solutions involve three kinds Bessel functions 
and since these are defined differently different books, perhaps worth 
stating the properties which define them here. Denoting differentiation 
with respect the argument dash, (assuming real 


and satisfy the differential equation 


The Bessel function arguments have been omitted some places below. 
Where the argument not specified implied that 
functions have argument 
functions have argument 
functions have argument pa. 


The easiest way proceed write down the solutions for and verify 
that they satisfy equations [18] and [19]. find: 
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Substitution from [108] [18] gives 


where the value for thermal neutrons the tank. Substitution from 
[109] [19] gives 


where the value for thermal neutrons the reflector. The coefficients 
and are fixed continuity current and density Assuming 
equal diffusion coefficients find: 
[113] 


can now evaluated from equation [17]. Omitting factors common 


therefore 
Io(qr) 
0 


Following the method section express the critical condition [106] 


from which obtain the differential relation (omitting Bessel function argu- 
ments before) 


Substituting from [21] and [22] for and gives 
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Replacing (equal thermal diffusion coefficients) and 
1)/T gives 


The first term predominates over the two other terms. Since 
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DEPENDENCE INTEGRATED DURATION METEOR 
ECHOES WAVELENGTH AND 


ABSTRACT 


The reflecting properties short-duration and long-duration echoes from 
meteor trails are examined. For short-duration echoes only, the observed relation 
between relative numbers and durations meteor echoes independent wave- 
length, the time day, the presence strong shower. Integrated duration 
times are determined from back-scatter experiments 9.22 m., 5.35 m., and 
2.83 m., and are found vary with the 3.5th power the wavelength, for either 
short- long-duration echoes. The integrated echo power depends the 6th 
power the wavelength. The effect changing the equipment sensitivity 
considered. The data from the back-scatter observations are used predict the 
integrated duration times the forward-scatter case, and particular, 
shown that the results the Cedar Sterling very-high-frequency 
experiment may explained meteoric reflections. 


INTRODUCTION 

Investigations Stanford University have shown 
flections from meteoric trails can support, intermittently, high frequency 
radio transmissions well beyond the line-of-sight, without the aid the usual 
layer reflections. The Stanford workers have developed theory 
this forward-scatter transmission, confirmed some experimental evidence, 
which indicates that the forward-scatter echo duration, can many 
times long the back-scatter duration, provided that the trail axis 
suitably oriented each case. the angle between the incident and 
reflected ray, they show that and this effect should greatly 
increase the percentage time during which forward transmission can 
supported compared back-scatter normal radar reflections. 

their analysis, the Stanford workers have assumed that each meteor 
trail long, thin, ionized cylinder, enduring for not more than second so, 
reflection highly dependent the orientation the meteor trail. 
The trail should tangent one family confocal ellipsoids, having the 
transmitting and receiving stations foci. the trail also lies plane 
containing the two stations, the amplitude the forward-scattered signal 
increased sec which additional factor favoring forward transmission. 

this paper shall use some the back-scatter echo data obtained 
Ottawa establish relation between the wavelength and the integrated echo 
duration time. From this, similar relation may computed (based some 
simplified assumptions) and used predict the performance the forward- 
scatter case. First, will pointed out that while the scattering from short- 
duration trails usually highly directional, the long-duration echoes tend 
lose this aspect sensitivity. 


received March 10, 1954. 
Contribution from the Radio and Electrical Engineering Division, National Research 


Laboratories, Ottawa, Canada. Issued N.R.C. No. 3305. 
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ASPECT SENSITIVITY 
have experimental data long-range km.) point-to-point 
communication via meteoric ionization. have observations forward- 
scatter over very short base line, obtained during couple hours the 
evening Aug. 11, 1949, when the Arnprior and Carleton Place radar stations 
were inadvertently operating nearly the same frequency, 


© 


205 25° 305 355 
AUG. II, 1949 


Fic. record meteor echoes Ottawa (A), Arnprior (B), and Carleton 
Place (C), Aug. 11, 1949. The forward-scatter echoes the Arnprior and Carleton Place 


records undulate slightly range. 


~ 
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which was also close the frequency the Ottawa radar. The transmitters 
Arnprior and Carleton Place each radiated kw. pulses, rate 120 per 
sec., locked the common 60-cycle mains supply. Consequently, when the 
phasing the two modulators happened suitable within narrow limits, 
the display each station would show not only the normal back-scatter echo 
but also quasi-stationary forward-scatter echo from the other station. The 
Ottawa transmitter was keyed nonsynchronous modulator about 
118 pulses per second, hence the echoes from Ottawa appeared the displays 
the other two stations vertical streaks about second apart. 

Fig. shows typical echo recorded all three stations 9.2 m., 
20" Aug. 11, 1949. The forward-scatter echoes can easily dis- 
tinguished from the back-scatter echoes the Arnprior and Carleton Place 
records the small sinusoidal wobble range, presumably caused 
slight phase variation the 60-cycle mains between the two stations. 

The interstation distances Ottawa—Arnprior km., Arnprior Carleton 
Place km., Carleton Place km. These distances are far too 
short enable draw any conclusions about the effect; one would 
not expect exceed more than for this reason alone. 
However, select meteor trails that have produced back-scatter echoes 
each station, then may expect the average significantly greater 
than the average This follows assume line-type scattering, with the 
geometrical consequence that the scattering point the trail for the forward 
transmissions will located between the two back-scatter points. the 
ionization density along typical meteor trail decreases either end 
the trail, the density the forward-scatter point will tend greater than 
the mean the densities the two back-scatter points. 

The Arnprior and Carleton Place transmitters were locked suitable phase 
for total time about two hours the interval 19" 20™ 22" 30™, Aug. 11, 
1949. The criterion echo selection was that the back-scatter echoes appeared 
all three stations, making sure the case very-short-duration echoes 
that the relative times initial appearances were consistent with the trail 
geometry. For the present purpose, only the Arnprior and Carleton Place data 
are used, but the Ottawa radar range was needed for triangulation and path 
determination selected meteors. Table summarizes the data for echoes with 
back-scatter durations less than one second either Arnprior Carleton 
Place, and for echoes with durations over four seconds. For each meteor the 
forward-scatter durations measured both Arnprior and Carleton Place 
were averaged, and similarly the mean the back-scatter durations observed 
both stations was used. 


TABLE 
RATIO FORWARD- BACK-SCATTER ECHO DURATION 
Echoes less than Echoes more than 
sec. sec. 


1.48 
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The number observations small because the stringent selection 
conditions but the results are nevertheless statistically significant. Without 
exception the mean was greater than the mean for each the short- 
duration echoes, while for the long-duration echoes the mean varied only 
few per cent, one way the other, from the mean each case. Since 
the effect negligible our experiment, the 50% increase over 
for short-duration echoes can reasonably explained qualitatively assum- 
ing that the short-duration echoes come from long, thin columns ionizations 
that are less highly ionized near the ends. 

The approximate equality and for the long-duration echoes could 
imply that the ionization density the back-scatter points the trail 
about the same that near the middle. This may not seem unreasonable 
first glance, one considers that large meteors may create trails 100 km. 
long, and that there could less variation the density per km. 
for selected km. segment than the case small meteor with over- 
all trail length km. However, was observed that for F-type echoes 
lasting few seconds the appearances the echoes each display were 
usually simultaneous (i.e. within 0.05 sec., which was the limit time reso- 
lution between stations), contrast the easily measurable delays observed 
the case the short echoes. The implication appears that all stations 
simultaneously receive echo from the same point the 
trail. This could either sudden distortion uniform trail (3), 
short segment denser ionization (7). Evidence from the simultaneous 
doppler station records confirms this impression; most the selected short 
echoes showed the first Fresnel zone while very few the longer echoes did. 

The Ahb echo types, example which appears Fig. show the 
simultaneous appearance the echoes the three stations. other 
echoes there may several seconds delay between the appearance the 
h-echo and the various e-characteristics, but separate triangulation the 
components shows that they arise from discrete points the original trail, 
and therefore that the scattering from these points must substantially non- 
directional. The corresponds closely the path the moving 
meteoroid. Every meteor echo occurring near may said have associ- 
ated head echo, caused diffraction scattering from the forming trail, both 
when approaching and receding from the observer, but only for the brighter 
echoes can appreciable change range detected. The rapid decrease 
signal strength either side agreement with the diffraction theory 
these cases. other cases the behavior the head echo suggests 
ionization the head the meteor. The diffraction echo has its 
maximum intensity the point, whereas many the observed head echoes 
‘do not behave this manner and are more consistent with the moving ball 
type target, with echo amplitude independent position with respect 
fo. 

The highly-directive properties the short-duration echoes and the non- 
directive characteristics the long-duration echoes were shown statistically 
previous paper (6) analyzing the short and long echo rates during the 
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Perseid shower the radiant rose. When the radiant elevation passed 60° the 
short-duration echo rates fell off rapidly because the ranges, corresponding 
the points, exceeded the maximum recorded range. the other hand, 
the number long-enduring echoes steadily increased, roughly the sine 
the radiant elevation angle, one would expect these echoes were inde- 
pendent trail aspect. 

sum up, seen that the short-duration echoes, say those less than two 
three seconds 9.22 m., behave the meteoric target were line 
scattering source, and hence have predictable directional characteristics. 
Echoes longer duration tend lose this aspect sensitivity and echoes more 
than sec. long may, the average, expected scatter energy all 
directions, and some may the first instant their appearance. 

There little the discussion far that new, the intention has been 
emphasize from previous knowledge the difference the reflecting properties 
long- and short-duration meteor echoes. earlier paper (7) postu- 
lated the presence denser ionization along the trail attempt 
account for some the long-enduring echoes. Greenhow (3) has criticized 
our hypothesis the basis that the blob theory were correct the echoes 
should detected short-wave radar before being seen long-wave 
radar, and experiments show that this does not occur, general. Greenhow 
postulates that the equipment sensitivities should adjusted the inverse 
third power the wavelength, which believe prejudging the issue, i.e., 
equivalent assuming that long-duration echoes have the same proper- 
ties short-duration echoes. However, accept this arbitrary sensitivity 
ratio does follow that smaller target area required produce the mini- 
mum detectable echo the shorter wavelength; hence, the blob expands 
should seen first the short-wave set. The flaw the argument the 
assumption that, given instant, the meteor target must present the same 
effective dimensions both radars. That is, the electron density across 
diameter the blob would have constant, with abrupt discontinuity 
the boundary. Owing radial diffusion this unrealistic picture; 
several better models have been suggested, and, for example, gaussian 
gradient density seems more appropriate. the critical density inversely 
proportional the square the wavelength, follows that the effective 
target area presented the long-wave set may considerably greater than 
for the short-wave set. The actual values depend the particular model 
chosen, but analysis several reasonable hypotheses indicates that the 
ratio the target areas can more than sufficient ensure that the long- 
wave echo stronger, despite the inferiority equipment sensitivity. 


BACK-SCATTER INTEGRATED ECHO DURATIONS 


earlier paper (6) presented some statistics back-scatter echo 
durations observed pairs simultaneously 9.22 and 5.35 m., and 
9.22 and 2.83 can use these data directly determine the integrated 
durations 5.35 and 2.83 m., but further reduction the 9.22 records 
will needed because many meteors detected this wavelength did not 
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appear the shorter wavelength displays. The discussion the previous 
section suggested that arbitrary value two three seconds might 
chosen dividing line for the two groups 9.22 echo durations. the 
average the durations 5.35 and 2.83 will only 1/3 and 1/10, 
respectively, the 9.22 durations for the same meteors. Therefore, shall 
divide the 9.22/5.35 data Ref. into echoes less than sec. and echoes 
sec. and over, and the 9.22/2.83 data Ref. into echoes less than 
0.4 sec. and those greater. 

define the total echo duration divided the total observing time; 
may expressed fraction (which may exceed unity) seconds per 
echoes duration and will refer the back-scatter and forward- 
scatter cases respectively. Letters brackets following the symbols will 
used where desirable distinguish between the echo groups, follows: 
(S) short echoes, (L) long echoes, (A) all echoes. will used for the frac- 
tional time during which least one echo present. will general less 
than and will never exceed unity. 

Table here summarizes all data from Table III, Ref. Table the 
data for Film No.’s and only, from Table Ref. Table all data from 
Table VI, Ref. and Table all data from Table II, Ref. The normalized 
hourly durations are the observed durations multiplied correction factor 
equate the system parameters the 5.35 and 2.83 sets that the 
standard equipment. This factor will derived later the section 
the variation with system sensitivity. here, allowance 
made for the fact that the increase duration with given increase sensi- 
greater for the shorter echoes (5), and slightly different factors are 
used for the short- and the long-duration groups. 


TABLE 
5.35 PERIOD. MEAN TIME 1850 E.S.T., 1950 


Observed (sec./hr.) 5.9 
TABLE III 


5.35 PERIOD. MEAN TIME 0430 E.S.T., 1950 


Observed (sec./hr.) 8.9 310 319 
TABLE 


5.35 PERIOD. MEAN TIME 2330 E.S.T., 1950 


Observed (sec./hr.) 8.1 120 128 
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TABLE 
2.83 AVERAGE. 1950 


Observed (sec./hr.) 0.22 0.51 0.73 
Normalized (sec./hr.) 2.8 8.4 


Qu 


Counts were made the long and short echo durations using 


Film No. the minimum period (Ref. and Film No. the maximum 
period the 9.22/5.35 observations, and these are listed Tables and 
VII. Here, significantly less than Dg, and the values may de- 
termined directly subtracting the measured overlaps. The overlap most 
cases was the order but during the maximum period amounted 
for the long-duration echoes. 


TABLE 
9.22 PERIOD. MEAN TIME 1850 E.S.T., 1950 
Ta < 3 sec. TR > 3 sec. All Tz 
Observed (sec./hr.) 129 120 249 
Observed 128 117 243 
TABLE VII 


9.22 PERIOD. MEAN TIME 0430 E.S.T., 1950 


Observed (sec./hr.) 227 1358 
Observed (sec./hr.) 223 1160 


Except for Table VI, the tables show that greater than and 
the ratio increases the shower radiant rises. This depends some extent, 
course, the selected dividing lines between the groups but not very 
sensitive the choice these lines, within reasonable limits. 

Film No. (Table V1) the 9.22 part the record showed meteor 
echo lasting min. which did not appear all 5.35 This the longest 
duration that have noted our records; was overlooked previous 
paper long-duration echoes (8). This exceptionally long echo has been 
deleted from the data Table its inclusion would have seriously upset 
the statistics. The echo extended from 200 km. 260 km. for much its life 
and hence will have obscured some normal echoes. The masking would more 
severe for the long-duration group than for the short-duration echoes, the 
latter echoes tended occur short ranges this period. 


determine the wavelength dependence the experiments the 
previous section should ideally have been done with all three radars simul- 
taneously. Or, lacking this, the and the 9.22/2.83 records 
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should each have covered shower period and nonshower period. Regretta- 
bly, the 9.22/5.35 data are available only for the period the Perseid 
shower, and the 9.22/2.83 data were observed over many days which 
strong showers occurred. Fortunately, the mean all the 9.22/2.83 
records may taken represent the average performance, around the clock 
during nonshower conditions, and period when similar average conditions 
existed may selected from the 9.22/5.35 information. 24-hr. basis 
the average nonshower echo rate the 9.22 radar has been established 
about 400-500 echoes per hour. plot hourly echo rates versus time 
showed that the 9.22 over-all rates were about 450 echoes per hour the 
neighborhood 2100 E.S.T. during the August period. Interpolation the 
data the tables above was then deduce values for 9.22 
and 5.35 2100 E.S.T. These values, together with the 2.83 normalized 
figures, will adopted averages for mean standard nonshower con- 
ditions; are given Table 


TABLE VIII 
ADOPTED AVERAGE FOR DIFFERENT WAVELENGTHS 


9.22 150 400 550 
5.35 100 118 
2.83 2.8 5.6 


The data Table VIII are presented Fig. plotting the common 
logarithm expressed sec./hr., against the logarithm the wavelength 
For convenience, scales and percentage time have been 


MEGACYCLES PER SECOND 
300 200 180 120 10090 80 70 60 50 40 30 


LOG Dg (IN Sec /HR) 


LOG > (IN METRES) 


Fic. Observed values log log Dg(L), and log for three wavelengths, 
9.22 m., 5.35 m., and 2.83 Back-scatter case under average conditions. 
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added Fig. The three straight lines have been fitted least squares, and 
have the following equations: 


The exponent depends here our selection 9.22 and 5.35 
observing periods, but not very sensitive this choice. For example, the 
range 400 700 sec./hr. for 9.22 echoes yields the spread 
may therefore conclude that the fractional time during 
which back-scatter reflections are obtainable proportional exponent 
the wavelength that greater than the 3rd power but not likely exceed 
the 4th power. 

For convenience our later discussion shall adopt 3.5 the power 
for all three lines Fig. The proportionality factors Eqn. will also 
require some readjustment, which will based the approximate values 
read from the graphs The adjusted set equations is: 


The relative values the proportional factors for and depend 
the arbitrary choice the dividing lines for the short- and long-duration 
groups. Using different dividing lines the range 9.22 m., 
higher than 

not have data for wavelengths longer than 9.22 m., but the experience 
the Stanford workers suggests that the lines Fig. might extrapolated 
approaches and exceeds unity will approach unity. elementary analysis, 
based probability theory, suggests that the meteors occur random, 
that the fainter visual meteors tend occur more frequently groups, 
rather than being distributed entirely random (9). This has been confirmed 
for radar meteors some unpublished analysis progress Harvard, using 
our radar records. 


METEOR NUMBERS AND ECHO DURATION 


Fig. have plotted the logarithm relative numbers echoes (not 
corrected hourly rates) against the mean logarithm the duration, for 5.35 
echoes. The circled points are taken from Table III, Ref. and represent 
period minimum background and shower activity, i.e., 1630-2110 E.S.T., 
during the Perseid shower. The crossed points are taken from Table Ref. 
and represent period high activity for all background meteors and the 
long-duration shower meteors (0210-0650 E.S.T.). The solid straight line has 
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been drawn through the circled points, with measured slope —0.74. 
The dashed straight line was fitted the short-duration points only the 
crossed points; its slope —0.82. The long-duration groups the maxi- 
mum period deviate considerably from the dashed straight line because there 
were large number long Perseid echoes this period but relatively few 
short-duration Perseids. 


LOG N 


LOG Ts 


Fic. between the logarithm the relative numbers echoes and the logarithm 
their durations, observed 5.35 The circled points and solid line represent the minimum 
period, 1630-2110 E.S.T., the Perseid shower, and the crossed points and dashed line 
represent the maximum period, E.S.T., Aug. 11-13, 1950. 


Fig. similar relations between relative numbers and durations are 
shown, using the 9.22 data, which are summarized Tables and VII 
this paper, but not tabulated detail. The circled points refer the mini- 
mum period, 1630-2115 E.S.T., and the crossed points the maximum 
period, 0207-0657 E.S.T. Two straight lines have been drawn through the 


LOG N 


~2 ° ' 2 
LOG Ty LOG Ts 


Fic. between the logarithm the relative numbers echoes and the logarithm 
their durations, observed 9.22 The circled points and solid line represent the minimum 
period, 1630-2115 E.S.T., the Perseid shower, and the crossed points and dashed line 
represent the maximum period, 0207-0657 E.S.T., Aug. 12-13, 1950. 

Fic. Relation between the logarithm the relative numbers echoes and the logarithm 


their durations, observed 9.22 The observations were made during two nonshower 
periods, 2230-0130 E.S.T., June 14/15 and 15/16, 1950. 
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short-duration groups for each period, and again seen that the long- 
duration groups deviate from the respective lines. The slope the solid line 
(minimum period) —0.72, and that the dashed line (maximum 
period) ism 

Fig. shows the log N/log relation for average nonshower period. 
The data were taken from the 9.22 records for June 14/15 and June 15/16, 
1950, 2230-0130 E.S.T. The slope the line drawn through the short-duration 
echoes —0.72. The long-duration echoes fall below the line this case. 
appears that during periods little activity the slope the log versus 
log graph for long echoes becomes more negative than the slope the short- 
duration graph, and during periods peak activity the slope becomes more 
positive. The long-duration echoes are the ones more usually correlated with 
visual meteors. Hence, anomalous conclusions may drawn, unless care 
exercised, one were attempt compare these log N/log relations 
directly with Millman’s log graph (10), combined with the log 
relation given previous paper (4), where the visual magnitude. 

seen that the slope the log versus log relation practically 
independent echo rates and wavelength, confine our attention the 
short-duration echoes. shall adopt slope —0.75 representative. 
had happened that —1, then all groups would have had the same 
duration sum, analogous the relation constant mass per magnitude 
interval given Watson (13). 

have used scale for the-duration groups (or nearly for the 
very short echo groups), From above, where 
—0.75. the over-all number meteors from the largest observed 
down T,, then 


n p=0 


The integrated echo duration 
n n 


VARIATION WITH SYSTEM SENSITIVITY 


Let the sensitivity our standard system varied changing the trans- 
mitter power the receiver bandwidth, keeping the wavelength constant. 
define the sensitivity factor the logarithm the power sensitivity 
ratio. has been found that the over-all echo rates vary the square root 
the sensitivity ratio (4), log x/2. From above, doubled, 
means four-fold increase power for example, follows that each will 
also doubled. Furthermore, every echo that was recorded before will 
increased slightly duration (5), the increase being greater for the short- 
duration echoes. Both effects will increase but the increase numbers 
more important. 

our work arbitrary lower limit resolution 0.05 sec. has been 
set. doubt there are myriads echoes shorter than this the records 


oo 
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but they tend indistinguishable from noise spots, which are the order 
0.02 sec. This simply mechanical limitation the recording system; 
the wide bandwidth 250 will inherently allow very much better 
resolution, but are limited the spot size the cathode ray tube and 
the 120 c.p.s. sampling frequency. Our doppler records, for example, although 
using lower power and narrower bandwidth, have resolution millisecond 
better. For purposes this discussion minimum 0.05 sec. for any 
sensitivity factor will assumed. This means that the sensitivity in- 
creased, more echoes will appear the group sec., and those 
that were there will increase duration, with many them now falling within 
the next group. the dividing line for short- and long-duration echoes 
our standard system some the echoes the group below will pass 
over above These are still echoes from directional trails and should 
included D(S). Because the fixed lower limit the range summation 
for D(S) will thus increase decrease with similar increase numbers 
and durations will observed for the long echoes, and the range summation 
will again vary slightly with 

For our standard system have computed 
experimental data for each interval log 0.301 the range —1.20< 
log 2.10. The summation for D(S) extended over the first five groups 
log 0.30, corresponding roughly our selected and the summation 
D(L) included the remainder. The sensitivity ratio was then increased 
factor 100, for example, which increased the numbers each group 
and added new group D(S) because the mean had been increased 
factor The number groups the D(L) summation remained the same 
creation new group beyond log 2.10. 

Fig. shows the relative values and computed this 
manner and plotted against The graphs are not quite straight lines but mean 


LOG RELATIVE UNITS) 
Nn 


LOG SENSITIVITY RATIO 


Fic. Computed relative values log and log plotted against the logarithm 
the sensitivity ratio. 
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straight lines can fitted either curve and each has the slope 0.56. Hence, 
have 
[3] log (System 0.56 log (System I), 


where the log the sensitivity ratio System System Inclusion 
the effect increased duration has thus raised the exponent the sensi- 
tivity factor from 0.50 0.56. 

Our standard system has peak power 200 kw. (average power kw.), 
antenna gain 3.5, minimum detectable signal watts receiver 
bandwidth 250 and assumed use any wavelength between 
and bandwidth kc./sec. were acceptable, for example, the peak 
power could reduced 0.8 kw. with the same sensitivity, though with some 
depreciation the resolution, the rate transmission intelligence. 
However, has been assumed far that the echo reflections were usable 
down receiver noise level, which undesirable for practical purposes. 
signals db. above noise are regarded the lower working limit, one could 
increase the transmitter power factor 10, reject all signals below db., 
and still have the same Or, alternatively one could leave the system para- 
meters fixed and use only echoes rising db. above noise. Here and 
the new 0.28 times the original Referring Fig. this would 
mean that the usable back-scatter transmission times would divided 
anywhere the chart. 


RELATION BETWEEN INTEGRATED ECHO POWER AND WAVELENGTH 
From Equation have for any two selected values and 
From Equation may write, 
where and are, for example, the two transmitter powers given 
wavelength that will yield the same and that would obtained with 
equipments identical sensitivity operating wavelengths and 
given Equation 
From Equations [4] and [5] may eliminate and hence, 
Thus, order maintain selected value the sensitivity must vary 
the sixth power the wavelength. the system sensitivity held constant 


then Equation [6] also indicates that the integrated echo power will vary 
the sixth power the wavelength. 


THE FORWARD-SCATTER CASE 


The results the back-scatter experiments may employed suggest 
the probable performance when the transmitter and receiver are separated 
the order 1000 km. Eshleman and Manning (2) have developed theory 
the probability detecting randomly oriented trails over such path. Meteor 


= 
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trails are not distributed randomly, certainly not their space orientations, 
and possibly not time. The distributions the known meteor radiants are 
very complex, and, will seen later, the meteors contributing the 
forward-scatter case may several magnitudes smaller than the meteors that 
have been studied hitherto and their distributions are unknown. therefore 
hesitate attempt extend the Stanford analysis more general case that 
might better representation actual conditions. fact, for the purposes 
the following discussion shall greatly simplify the working assumptions, 
and avoid probability theory and meteor distributions much possible. 

Following the Stanford experience shall simply assume that the observed 
may multiplied factor greater than unity, based the 
effect mentioned earlier, order yield computed Villard (11) 
found for reflection the mid-point 1175 km. path that 23. 
The average ratio for all distinguishable echoes was 4.5, which included many 
cases where the reflection point was near one station and was therefore 
much smaller. shall assume that directional antennas will used, 
illuminate only region the vicinity the mid-point. Thus, the average, 
the ratio should intermediate between the above two observational 
values. shall adopt representative multiplying factors, f(S) 10, and 
15. The amplitude factor, proportional sec will also operate favor 
increased but shall amit this for the moment. will not 
increased the effect much but the forward-scatter dura- 
tions should certainly increased some factor because many the long- 
duration trails will have retained some directional characteristics. Arbitrarily, 
shall assume representative factors f(L) and 

Three computed graphs versus wavelength are shown Fig. 
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Fic. Computed graphs log function wavelength, for three arbitrary 
multiplying factors, f(5, 2), 4), and f(15, 6). Forward-scatter case under average 
conditions. 
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for for example, means that the calculated forward-scatter 
duration has been obtained multiplying and and 
then adding the two figures. The relations Eqn. have been used the 
basis. Because the assumptions involved, the graphs Fig. should 
regarded only representative the order magnitude expected 
the forward-scatter case, though believed that they are conservative 
underestimates. 

The graphs Fig. apply spaced-station system using the 
power and sensitivity our standard radar. The for system dif- 
ferent sensitivity and wavelength may found first entering Fig. 
the desired wavelength and then correcting the corresponding using 
Fig. the equivalent analytical formula, Eqn. adjust for the difference 
system sensitivity. frequency Mc./sec. for example, would 
expect, under average nonshower conditions with our standard system, that 
D,(A) 16, 24% the time, corresponding the assumed factors 
2), 4), and f(15, 6). During the early morning hours, with shower 
present, the values would roughly doubled, while the early evening hours 
they might halved. The presence strong shower might expected 
increase the time percentages further, when the radiant favorably located 
with respect the transmission path. 

The geometry the forward-scatter case differs from that the back- 
scatter experiments, and examination some the differences necessary 
insure that the above analysis justified. Let assume that the trans- 
mitter and receiver are 1000 km. apart, each using similar directive antennas 
with patterns illuminating the large area centered above the 
mid-point the stations. The effective area the meteor region that typical 
20° 30° beam thus intersects about the same that intersected overhead 
our back-scatter radar with broad antenna pattern, when the recorded 
range limited 260 km. The relative numbers meteors given magni- 
tude will therefore comparable. The forward transmission path roughly 
three times the mean back-scatter path, and the attenuation with range will 
vary for the short echoes and for the long echoes; the forward 
loss will then 14-19 db. The antenna gain the 20° 30° beam (measured 
between half-power points) about 45, 16.5 db., compared the gain 
3.5, 5.4 db., used the back-scatter experiments. The net gain twice 
the difference these antenna gains expressed db.; this case the net 
gain will db., more than sufficient compensate for the range attenu- 
ation. Therefore, our estimates above may regarded reasonable 
and conservative. Admittedly, there considerable speculation throughout 
this analysis, especially the factors. There really satisfactory 
substitute for performing the actual experiments. 

practice, would probably better direct the antenna beams one 
side the other, both sides (split beam) the line joining the stations, 
Eshleman and Manning have recommended (2). Only meteor trails nearly 
parallel the earth’s surface and the line joining the stations will show the 
full benefit the and sec effects, and very few these occur. 
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directing the beams one both sides, greater number suitably-oriented 
trails will available from meteors with higher radiants. The long-duration 
echoes are not affected the same extent this reasoning. 
THE CEDAR RAPIDS-STERLING EXPERIMENT 

another example, consider the data forward-scatter experiment 
which the transmitter, 49.8 was located Cedar Rapids, lowa, 
and the receiver was Sterling, Virginia (1). The parameters this experi- 
ment, here designated System are compared with our standard character- 
istics evaluated for the forward-scatter case, System and the net gain 
loss shown Table IX. 


TABLE 
standard Cedar Rapids Sterling Gain db. 
Transmitter power 200 kw. (peak) kw. 
Path distance 1000 km. 1270 km. 
Receiver bandwidth 250 +42 


Net db. 


The receiver bandwidth for System has been calculated the geometric 
mean the published values I.F. and output bandwidths, which were 
and 1/12 respectively. Here 3.0, and therefore for 
for System will roughly 385, 770, depending the choice 
97.9, 99.95, and for the respective choices. These figures apply 
signals measured down noise level, and suggest that highly improba- 
ble that the signal will fall below noise level for any significant time. 

For signal level db. above noise, which appears representative 
mean value the signal strength recorded System multiply the 
figures above yield percentages 108, 216, and 324 
respectively. The corresponding values are 66, 88.5, and 96.1%. These 
should regarded minimum figures for average conditions, because 
have ignored the contribution that will made System echoes less 
than 0.05 sec. duration. These may the order milliseconds and still 
pass through the I.F. amplifier; the output filter, course, integrates the 
power contributions all echoes less than few seconds duration. The over-all 
-echo rate System will times the rate System over 200 echoes 
per min. duration 0.05 sec. and longer. With time constant sec. 
the output stage System not surprising that the signal should appear 
out and longer gaps should rare indeed. System will detect meteors 
least four magnitudes fainter than System possibly down the 14th 
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15th visual magnitude (4). The presence the strong visual meteor showers 
should not affect System much System because these showers 
appear contain relatively few very small meteors. 

The mean diurnal variation signal strength published for System 
the same general form the typical diurnal meteor rate curve obtained 
with our standard back-scatter system, though there are some discrepancies. 
These may due the hitherto unknown distribution meteors the 
range magnitude. far are aware none the conventional 
high-resolution meteor systems has been able study meteors this magni- 
tude. Naturally, the individual identity the meteor echo entirely lost 
System except for its average power contribution. may that incorpo- 
rating filters with long time constants the present meteor systems will add 
our knowledge these small meteors, though the observing techniques and 
the analysis will have modified considerably. 

have shown for System that should very close 100%. Thus, 
continuous communication should obtainable, though low intelligence 
rate because the narrow bandwidth. Intermittent communication higher 
intelligence rates has often been demonstrated.* The authors the report 
the Cedar Rapids Sterling experiment considered that the received signal was 
due scattering from discontinuities the lower ionization, possibly 
due the effects turbulence. are suggesting here that their results might 
explained almost entirely the basis meteoric reflections. previous 
paper (8) reported the observation Mc./sec. weak scatter echo 
height about km., and that time suggested that the agency 
responsible for this continuous echo might the same the mechanism that 
supported the transmissions the Cedar experiment. 
this continuous echo were due meteors, the range meteoric masses 
must much smaller than the usual ranges background and shower 
meteors, and there can few meteors brighter than the 9th magnitude, say. 
Meteoric dust clouds could responsible. The case not proved and still 
feel that ionospheric discontinuities may offer stronger explanation for our 
weak scatter echo. the other hand, returning the forward-scatter case, 
think that meteoric reflections constitute the major contribution the 
observed signal, with scattering from ionospheric turbulence possible 
secondary mechanism. 
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EXPERIMENTAL STUDY BAND INTENSITIES THE 
FIRST POSITIVE SYSTEM 


VIBRATIONAL TRANSITION PROBABILITIES! 


ABSTRACT 
Integrated intensities bands the first positive system have been 
measured using recording infrared spectrometer. These data have been inter- 
preted vibrational transition probabilities which were compared with theo- 
retical values calculated under the assumption that the electronic transition 
moment independent internuclear separation. The comparison shows that 
the assumption not valid for this electronic transition. 
INTRODUCTION 
Accurate measurements intensities molecular spectra have large 
extent been neglected comparison with the establishment accurate 
wavelengths and energy level assignments. However, the study physical 
conditions molecular excitation existing spectroscopic sources both 
terrestrial and astrophysical origin, the need for accurate intensity measure- 
ments becomes evident. This may seen the dependence upon 
equation [1], 


2 


where intensity band recorded the spectrometer, 
population molecules excited vibrational level the 
source, 
energy separation between and levels, 
vibrational wave function the upper vibrational level 
vibrational wave function the lower vibrational level 
electronic transition moment, 
internuclear separation, 
aconstant depending upon geometry the optical system and 
the units employed. 


the vibrational transition probability the band con- 
cerned, and order that equation [1] may used, together with experimental 
measures band intensities, investigate, for example, the vibrational 
excitation set source, some prior knowledge the vibrational transition 
probability must obtained. The vibrational wave functions are readily 
accessible (3) for most molecular levels. The transition moment how- 
ever, presents more severe difficulties very formidable task calculate 
for each electronic transition interest. has been conventional assume 
that very slowly varying function (5) and thus that for each the 
bands 


Vv 
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where the average this slowly varying function over the contributing 
range Accurately computed values squares the overlap integrals 
are available for number important band systems (4, 7). 

importance therefore make experimental investigation into the 
validity the assumption constancy order determine whether 
the calculated values are sufficient knowledge the vibrational 
transition probabilities, whether any allowance for the variation with 
will have made. This paper describes such experimental investigation 
which experimental values have been used obtain relative values 
for bands the first positive system using equation 
comparison between these values and those the overlap integral 
squares shows that, fact, some allowance for the dependence 
upon will have made, and that the simplifying assumption con- 
stancy too crude this case. The evaluation the dependence upon 
described the next paper (13). 

Until the present work, set reliable band-intensities for the first positive 
system nitrogen has not been available. Some workers 
(10, 11, 12) have reported peak intensities for large number bands while 
others (1) have made careful measurements the integrated intensities the 
bands the sequences only. the work described here, the 
relative intensities bands belonging the sequences —3, —2, 
with the overlap integral squares calculated Jarmain and Nicholls (8). 

EXPERIMENTAL 

The instrument used making the intensity measurements was Perkin- 
Elmer Model 12-C infrared recording spectrometer. was equipped with 
dense flint prism and lead sulphide cell. The dispersion this instrument 
the exit slit was approximately 150, 580, and 920 A/mm. wavelengths 
6000, 10,000, and 16,000 respectively. The lead sulphide cell was tested for 
linearity and was only used the region intensities where its response was 
found linear. The normal collimating system the spectrometer was 
supplemented quartz lens front the discharge tube, since occupied 
different position from the globar source normally used with the spectrometer. 

The spectrum the first positive system was excited water-cooled heavy- 
current capillary discharge tube the Pearse-Hunter (6) type, but somewhat 
larger than the original. The discharge ‘capillary’ had internal diameter and 
length 1.2 and cm., respectively. This ‘capillary’ was viewed end-on 
through quartz window the spectrometer. 

The discharge was maintained direct current generator capable 
supplying maximum current amp. 2000 The output voltage the 
generator was stabilized electronically, and this was necessary order 
maintain the light output from the tube constant for the whole period 
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recording the spectrum. After turning the discharge and adjusting pressure 
and voltage conditions obtain almost constant current, the intensity 
typical band (e.g. (0,0)) was recorded during the five-hour warm-up 
period necessary establish equilibrium conditions the tube. drift 
intensity few per cent was observed during this time. After this warm-up 
period, the variation the deflection shown the recorder was about 
+0.5%. will noted that this figure includes any variations the am- 
plifying and recording system well any fluctuations the intensity out- 
put the tube. was under these last conditions that the spectrum was 
recorded, with further adjustment tube conditions (1220 v., 0.9 amp., 
1.25 mm. Hg). These conditions were found optimum for intense 
radiation the first positive system. 

The prepurified nitrogen used these experiments was obtained from the 
Mattheson Company, New Jersey, and was pumped continuously through the 
tube. The recording the spectrum was interrupted number wave- 
lengths ascertain that the intensity the (0,0) band, which was used 
‘internal monitor’ intensity variations, had not fluctuated during the 
recording. The slit width used was 0.013 mm. except for the very weak sequence 
5370-5630 for which the slit width was increased 0.025 mm. 

Intensity calibration measurements the optical system and detector com- 
bined were made recording the spectrum General Electric standard 
lamp with tungsten filament and quartz window, placed that the center 
its filament was the position previously occupied the near end the 
discharge capillary. This lamp had been calibrated the makers. The spec- 
trum mercury arc placed the same position was also recorded order 
provide data for wavelength calibration and dispersion correction for the 
standard lamp spectrum. The data obtained from these two spectra were used 
together with the tabulated values the emissivity tungsten interpolated 
from the results collected Forsythe and Adams (2) give values the 
relative efficiency the whole optical system approximately 257 intervals 
corresponding the wavelengths which the radiation function was tabu- 
lated (9). Values the relative efficiency versus wavelength were plotted 
large graph, and the relative efficiency for any wavelength could read from 
smooth curve which was drawn through the plotted points. was noted 
that there was practically scatter the values the relative efficiency. 
However, this was only partial indication the accuracy with which the 
values had been determined. The complete curve would have been tilted the 
temperature used determining the values Planck’s radiation function 
were not the true temperature the lamp filament. was found that the ratio 
radiant energy per unit wavelength interval 5400 that 17200 
(approximate range intensity measurements) would error approxi- 
mately 0.4% for each degree the temperature might have been error. The 
precision the manufacturer’s calibration the standard lamp not known 
but expected that not error more than few degrees. 
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REDUCTION DATA 


The profiles the bands observed were not ideal for immediate interpreta- 
tion terms relative intensities because the complicated blended rota- 
tional structure each band, and also because the partial overlap bands 
caused both development rotational structure and the dispersion 
the spectrometer employed. The following procedure ‘synthetic’ separation 
the band areas was therefore employed. 

For most bands, the sharply falling leading edge each band was continued 
straight line down zero intensity and the shape the tail the adjacent 
band was then found subtraction ordinates. For few weak bands each 
which appeared blended the tail very strong band, was necessary 
reverse the procedure and continue the tail the strong band down 
zero intensity and subtract ordinates. Although this procedure open 
some criticism, the simplest objective one which can used when the 
bands are not completely separated. The area under each band separated 
was measured with polar planimeter and then corrected for the relative 
efficiency the system from the calibration curve described above. The rela- 
tive intensities obtained this way are shown scaled (0,0) 1000. 


TABLE 
RELATIVE INTENSITIES OF THE BANDS OF THE FIRST POSITIVE SYSTEM 
” 


INTERPRETATION RESULTS 
From equation [1] the ‘vibrational transition probability’ may written: 
Poo’ = Ty or? /GNy 


where constant depending upon the geometry the system and upon the 
units and the wavelength corresponding the center the 
band half intensity. Although usual calculation this sort use 
the wavelength the band head, the center the band half intensity seemed 
more reasonable objective choice here view the fact that the bands 
not exhibit marked heads and band shape changes from band band. 
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The wavelengths which were used are shown Table II, and Table III lists 
that since the values for each (v’ const.) still contain 


TABLE 
WAVELENGTHS OF THE BANDS OF THE FIRST POSITIVE SYSTEM IN A 
6,824 7,569 8,504 9,594 11,110 12,844 15,370 
6,098 6,661 7,338 9,147 
5,496 5,931 6,432 
5,424 5,829 6,291 
5,352 5,730 
TABLE III 
Pore No FOR FIRST POSITIVE SYSTEM 
| | 
663 434 570 736 
419 173 150 351 477 
26.4 235 197 133 227 479 
8.69 85.2 147 108 
3.70 15.8 
| 


comparison between transition probabilities from theory and experiment can 
only made along v’’-progression, for each band which will the 


same. 


For each the values given Table III have been 
rescaled and are compared Table with the theoretical squares overlap 
integrals previously referred (8). For ease comparison, the scaling has 
been such that entries both theoretical and experimental data for the first 
band each progression are equal. 

order allow comparison easily made between these two sets 
data, they are presented block diagram form Fig. will noted from 
this diagram that while there some agreement between the relative sizes 
the experimental and theoretical values the transition probabilities for 
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COMPARISON OF THEORETICAL AND EXPERIMENTAL TRANSITION PROBABILITIES FOR BANDS 


v’ | 0 1 
0 .340 .323 
(.340) | (.282) 
1 .406 .002 
(.406) 
2 . 200 212 
(.200) | (.388) 
3 .301 
(.050) (.445) 
4 .006 .134 
(.134) 
(.027) 
6 -003 
(.003) 
| 
11 | 
12 | | 


.103 
. 266) 


to 
ore 


(.521) 
-210 


(. 265) 


059 
(.036) 


.008 


(.008) 


OF THE FIRST POSITIVE SYSTEM 
148 .089 .042 .018 .009 .004 .002 
.002 .074 126 119 .078 .043 .021 .008 
031 .009 .069 .097 .089 .076 .041 
(.088) (.377) 
.046 .044 104 .037 .000 .023 .057 
.259 .083 .107 .003 .084 .074 .014 .002 
012 .116 .267 .030 .107 
! | (.116) | (.495) | | 


11} 


.073 


-000 


12 


-008 


.056 


Note: The experimental values are given parentheses scaled such that the value for the first 
band each equal the theoretical value for the same band. 


positive system. 


Comparison theoretical and experimental transition probabilities for the first 


| .190 
| (.213) 
| .005 | .002 
| .021| .012 
| .055 | .036 
.073 | .063 
| .028 .051 
| .003 | .006 
-061 | .017 
-033 | .063 
.014| 
| 
1 1 
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neighboring bands nevertheless the experimental and 
theoretical values tend diverge along the progression. Passing out along 
progression higher quantum numbers entails the use increasingly 


and Yor noticed for increasing range integration suggests that 
the simple assumption being sensibly independent not justified. 

The data presented here will used the following paper find the depend- 
ence the electronic transition moment internuclear separation. 


The research reported this paper has been sponsored the Air Force 
Cambridge Research Centre under Contract AF-19(122)-470. 
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EXPERIMENTAL STUDY BAND INTENSITIES THE 
FIRST POSITIVE SYSTEM 


II. THE TRANSITION MOMENT! 


ABSTRACT 


Experimental band intensities the first positive system 
nitrogen are used together with theoretical transition probabilities determine 
the dependence the electronic transition moment the internuclear separa- 
tion. Over the range 1.17 1.59 internuclear separation, this depend- 


INTRODUCTION 


part hereafter referred (3), measure agreement was reported 
between the relative magnitudes theoretical (overlaps squared) and experi- 
mental vibrational transition probabilities for neighboring bands 
progression. However, there was decided trend experimental values away 
from the theoretical values along This divergence suggested 
the invalidity the conventional simplifying assumption made computing 
the vibrational transition probabilities, namely that the electronic transition 
R,(r) for the system was roughly independent internuclear separa- 
tion this assumption, the vibrational transition probabilities had been 
computed the squares the overlap integrals the corresponding wave 

this present paper, experimental intensities and the theoretical overlap 
integrals are used conjointly obtain the dependence R,(r) for the 
first positive system Such information very difficult obtain 
theoretical means, since would require detailed knowledge the electronic 
wave functions the states involved. With knowledge the dependence 
the previous overlap squares may modified give array 
which represents more closely the values the vibrational transition prob- 
abilities. These are indicative the distribution the ‘oscillator strength’ 
the electronic transition through the bands the system. 


METHOD AND RESULTS 


The method used here for determining the form the electronic transition 
moment has been suggested and fully described Fraser (1). shows that 
the following may written for the ratio the values the electronic tran- 
sition moment for any two bands 

[1] R. (Fev) / Males 

equation [1] and are characteristic wavelengths for the respective 
bands while and are respectively the integrated intensities the 
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same bands. The subscript refers one particular band v’’-progression. 
Which kind average internuclear separation for the tran- 
sition, defined by: 

(v’, r, v’’) 


where 


Equation [1] enables experimental determination made the vari- 
ation the ratio along using the experimental 
results and the previously reported theoretical values (2). These 
ratios for each observed v’’-progression the system are shown 


TABLE 


FOR BANDS THE FIRST POSITIVE SYSTEM 
(the scale different for each 


| 
1.00 1.61 1.41 1.74 
1.00 1.35 1.19 10.3* 2.09 1.86 
1.04 1.00 1.26 2.19* 1.45 
| | 


order interpret the variation function the value appropriate 
each band the system must found. The were calculated the 
simple method suggested Jarmain (1) rather than through the definition 
equation [1]. The for large number bands are given Table 
and thus value may associated with each ratio 


TABLE 
BANDS THE FIRST POSITIVE SYSTEM 


eo | | | | | | | | 
| | | | | | | | | | 
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Before discussing how the ratios Table are brought the same scale and 
thus linked together over the whole range covered the the ob- 
served bands, brief mention will made the reasons for rejecting some the 
points which this stage were considered unreliable. For each 
graph versus was plotted, and the data for points which 
lay far from smooth curve representative the trend the other points were 
discarded after careful examination. These bands are marked with asterisk 
Table The discarded bands had low overlap integral all cases, and 
was thought that the percentage error these might relatively high. Also 
there may have been some overlapping the sequences the spectrum. The 
resultant overlapping band one sequence with weak band from another 
sequence would probably not noticed some cases change profile 
the low dispersion the spectrometer used. Thus the intensity taken due 
one band the sequence might actually due both and 
weak band far out the sequence This would introduce large 
error only the case weak bands. 

The remaining points each graph versus were con- 
nected straight lines. The relative values from two 
progressions were brought the same scale normalizing the areas under 
the lines the common interval This process, which demonstrated 
was used set the data for all progressions the same scale the 


15 

0.5 

| 

130 1.35 1.40 1.45 1.50 1.55 


Values progression have been brought the same scale the 


| Fi 
| 
0.5 
| 
| | 
0 'H 
1.35 1.40 1.45 1.50 
1.5 
oo 
1.0 
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progression. These scaled values are given Table 
They represent the variation R,(r) with for the band system. While 
the values the Table III are all expressed the same scale, 


TABLE III 
FOR BANDS THE FIRST POSITIVE SYSTEM 

0.653 1.05 0.920 1.14 

0.424 0.585 0.776 


this scale still arbitrary. Absolute values could obtained only through 
theoretical calculations lifetime measurements. 


DISCUSSION 


attempt represent the information given Tables and III 
functional form, the values were assumed exact, and was found 
that 
was good representation the dependence over the range 
covered the bands the system that were used (1.17 1.59 A). 
The constants equation [3] were found from least squares 
log versus Because the arbitrary scale R,(r) significance can 
attached the constant factor equation [3]. The curve and the data 
which represents are shown Fig. Following standard procedure, the 
probable error the ratio two values 


exp (3.61 3.02r) 


1.4 


rin Angstroms 


1.3 


Fic. versus for the first positive system nitrogen. 
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Using equations and [4], the following values and limits are obtained: 


will noted that these limits have been determined the scatter the 
data. However, was pointed out any inaccuracy the manufacturer’s 
calibration the standard lamp used the work would further widen the 
ranges error given above. 

Thus expression has been found represent the shape the electronic 
transition moment for the first positive system nitrogen. This may used 
calculating more accurate values vibrational transition probabilities. 
However, must noted that can only used for bands whose falls 
the range covered the values used determining equation [3], 
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FUNCTIONAL POTENTIAL 


ABSTRACT 
covariant two-particle wave equation the following form investigated: 


both indices. the internal wave function depending only the relative 
coordinates. This equation provides covariant model which exhibits nonlocal 
interactions and can studied relatively simple methods. The investigation 
primarily methodological. The physical model similar the pion 
and like it, qualitative and not based fundamental theory. 


INTRODUCTION 

has been suggested Fermi and Yang (6) that the pion consists 
nucleon bound antinucleon.* Since such model does not fit very well 
into the framework current field theory, any elaborate investigation 
terms this theory would probably not promising; but course other 
theories the situation may different. this paper such alternative theory 
not considered, and the existence the composite pion simply assumed 
without any attempt justification terms fundamental theory. One 
conceivable structure then investigated means formally simple 
wave equation which covariant and has for bound state problems approxi- 
mately the same physical content the picture. 

This equation 


the eight coordinates and transforms spinor both indices; 
responsible for the interaction between the two particles and formally resembles 
potential well. order give simple discussion Fermi and Yang assumed 
square well interaction; equation [1] may regarded having been ob- 
tained from their picture replacement the Dirac hamiltonian plus 
square well covariant operator. 


received April 1954. 
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The preparation this paper was sponsored (in part) the Office Naval Research, United 
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*There were earlier suggestions the same nature Broglie (2) and Rosen (11). The work 
reported here was finished about the same time three other studies the same model: Broglie 
and Tonnelat (4); Broglie (3); and Moseley and Rosen (10). the present time the combination 
equation [1] would appear more interesting than but the same method may 
used the two cases. Unfortunately the multiplicative combination has not been studied. 
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chosen simple invariant which becomes important when the two 
particles are close For example and has the range 
Mc, then the effective potential well will have comparable but somewhat 
smaller range. The constant has the dimensions Fermi constant and 

Although [1] covariant, may deviate from relativity, since 
depending its form—may permit interactions outside the light cone; 
thus equation [1] general enough exhibit nonlocal interactions. This 
action-at-a-distance goes beyond the correlations already present quantal 
descriptions (5),* especially does not depend explicitly coordinate 
separation. the following, since shall concerned only with bound 
state problems, may ignore any explicit dependence coordinate 
separation without getting unphysical results. 

TRANSFORMATION PROPERTIES 

require that transform like the probability amplitude for 
and denote spinors belonging positive energy neutron and negative 
energy proton respectively. Let where charge conjugation 
operator; then transforms like where the wave function 
positive energy antiproton. characterize the (negative) pion 

Under the Lorentz rotation 


Under space inversion 


where for example, depending upon inversion properties assigned 
Take both and type (Yang and Tiomno (12)) that 


follows: 
0 1 2 1 


clear that equation [1] not proper quantal equation since the probability amplitude does 
not even satisfy linear equation. There some formal resemblance between the model studied here 
and the work Bohm (1) that the wave function plays dual role, being responsible for forces 
well for probabilities; but the similarity stops here. 
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where totally antisymmetric and absolute value unity. 
Write 
‘ 1) 2 


where stands for all the subscripts 
Under the Lorentz rotation 


Since have taken follows from [7c] that the p-compo- 


The treatment the adjoint functions and (transforming like 
and respectively) follows similarly the 
definitions 


[8d] yt= 


THE EQUATIONS MOTION 


The equations motion may derived from the Lagrangian 


where 
[9a] 


[96] Lo = (3) + Pu” an 


[9c] 
and the interaction term. terms this Lagrangian 

where the number times appears thus 


and had been assigned opposite type would have been the p-component 
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The equations motion 
become 


‘ 


The composite system separated into and 


one the usual transformation 
(2) (2) 


Xu (3) (x, Xp h Xp = Xp — 


where are the external (internal) momentum operators; and 
are the corresponding position coordinates. the case the free pion one 
may further write 


o 2), 1 o i 


and the equations motion become 

where 


will now explicitly assumed that the interaction depends only the 
internal wave function, therefore one may put The equations 
motion can then written the following form: 


The notation used here does not distinguish between covariant and contra- 
variant components. 
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SPECIAL SOLUTIONS THE PROPER SYSTEM THE PION 
The probability amplitudes determined the set [16] depend upon the 
times both particles general. There are, however, special solutions which 
not depend upon the relative time have been able investigate 
only these special solutions. These solutions may defined the covariant 
condition 


the proper system, for which im, [17] becomes 
[18a] 


that there dependence the relative time. For these solutions one 
may write 


where the right-hand side the times particle and particle have been 
set equal. 

may noted that [17] satisfied the limit interaction, for 
then 


Although there reason require [17] general, this condition may 
used further limit the model; that will the procedure here. 

The solutions have been able obtain are also characterized simple 

0 1 3 

angular dependence the proper frame for the pion, namely: Xo, 
are either spherically symmetric vanish, while the remaining functions 
have p-orbital dependence vanish. Such angular dependence corresponds 
composite particle with zero spin (which correct for the physical 


pion). 
These solutions may written the following special form: 
do = A(r), = H(r) (r:/r), BIN = Aa(r), 


The equations motion the proper system are then: 
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imaginary the center mass system. These equations split into two sets 
involving (a, and (H, A), respectively, coupled only through the 
interaction term, simple solution the complete set obtained 
put 
and solve the remaining equations for and with the following 
set coupled equations: 

THE INTERACTION TERM 
For simplicity the interaction term, was chosen the following form: 


[23] 


(When factorable, may written the familiar form: 
p 
the notation introduced for the center mass system one has 

Tz 4= — 4A? — 2? — 4A?, 


actions. 
SOLUTION THE EIGENVALUE PROBLEM 

The structure the composite particle determined the set [22]. 
Solutions are physically admissible only they are everywhere regular and 
quadratically integrable. This condition leads relation between the two 
parameters, the coupling constant, and the mass the compound system. 

convenient first write [24] dimensionless form. Put 


The dimensionless equations are then 


[25d] dx + vA + + On, 0, 


p 
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Eigensolutions exist only for particular relation between and The eigen- 
problem will now discussed method described more detail elsewhere 
(8)*. One may begin letting correspond time and investigating the 
motion the representative point the (It not necessary 
consider motion the since may eliminated the use 
the algebraic equation [25c].) then convenient consider the related 
obtained deleting terms which contain explicitly. 
This system 


9 > — = 


This set has the denoted and given 


For the conservative system [26] 


but for the complete set [25] 


The representative point for [26] thus moves the curves constant, 
while the corresponding point for [25] crosses these contours—in the direction 
larger smaller values depending the sign Note that 
order draw the one must eliminate from [27] 
[25c]. Note also that this plane eigensolution defined curve 
which starts and terminates the origin without ever going 


infinity. 
now consider the vector interaction detail. Take 


The integral which will called the contour function, 


and 


*The physical problem there completely unrelated this paper. 


ay 
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Since the representative point moves toward higher values 
point moves outward instead toward the origin; therefore cannot describe 
eigensolution. fact corresponds toa repulsive interaction and for 
bound state solution require that Therefore take and put 


The differential equations [25] then become 


+E=0, 
addition one has 


may eliminated use the cubic equation but course one must 
choose the root belonging eigensolution. This choice may made 


| 
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constant contours the plane. (Shaded area denotes the region 


shown Fig. from which one concludes that solution terminating 
the origin (such eigensolution) must correspond the branch 
With this information one may then plot the contours the plane 
(Fig. 2). clear from this figure that representative point which starts 
anywhere the G-axis always moves toward higher until finally ap- 
proaches one the three points P_, the origin—these are the only 
possible end points. enters either the loops (shaded) will 
the corresponding maximum, either eigensolution must there- 
fore avoid the shaded area order reach the origin. The points P_, 
and the origin represent constant solutions the differential equations; 
linearizing the neighborhood these constant solutions, easy show 
that the representative point spirals into either the other hand the 
origin saddle point. 

With this qualitative information one may obtain the exact eigensolutions 
the following numerical method. One finds numerical integration that 
one solution starting (Go, ends and that another starting 
goes P_, and that neither curve goes below the G-axis. narrowing 
the interval one may then find the eigensolution any desired 
accuracy; and neither curve goes below the G-axis the corresponding eigen- 
solution will nodeless. moving the starting point the right one obtains 
eigensolutions with more nodes. 

Thus the starting value may take only discrete set values for 
given (On the other hand, one may hold fixed and vary then only 
particular values will correspond eigensolutions.) the present problem 
the mass the pion and experimentally fixed. addition fixed 
the probability normalization terms Holding fixed, one varies 
(or until eigenfunction obtained. The eigenproblem thus establishes 
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relation between and The procedure here assign its experimental 
value and then calculate and other quantities interest. 


NORMALIZATION AND RESULTS 
The following normalization condition was used: 


which simplifies for the solution obtained to: 
0 


With the parameter 0.08 the zero node solution was obtained with 
differential analyzer. These solutions are plotted Fig. where G(x) and 
E(x) are shown. The probability density shown Fig. 

phase plane analysis shows that for the value the 
parameter eigensolutions exist for the interactions Jo, and 


—E(x) 


G(x) 


0.5 


Fic. Plot G(x), E(x) obtained the differential analyzer. 
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1.0 a(x) 


0.5 Ke) LS 
x 


Plot p(x) —16F(x) G(x), using G(x) from Fig. and F(x) calculated from 
eq. 
Note: From the plot F(x) G(x) 1.20. 


CONCLUSIONS 

Although there theoretical justification for the physical model, still 
because its mathematical simplicity have used estimate the masses 
the neutral and charged particles and their lifetimes for and decay (9). 
(The initial wave function the two-nucleon system taken from the above 
calculation. The one-step process and the y-decay involves two 
steps; the calculation therefore requires one lower order perturbation theory 
than the conventional meson calculation. For example, y-decay the com- 
posite pion pictured just the same way the y-decay positronium, and 
does not require preliminary dissociation into heavy particle pair.) The 
calculated mass difference between the neutral and the charged particles, 
well the lifetime for correct, but the y-lifetime 10° times too 
short. These results, however, depend mainly the size the localized 
wave function, not its form, and hence tell nothing about details the 
model.* investigate these details one must reconsider the y-lifetime, and 
one must also apply the same theory the deuteron. pointed out 
Fermi and Yang the deep potential well the pion must over into high 
repulsive barrier the deuteron. That occurs only for certain choices the 
interaction the present model also. 

result its formal simplicity the fundamental equation [1] may 


discussion the composite pion not employing particular model has been given elsewhere 
(7). equation (6) that paper the inequality should reversed, but the form, even existence, 
universal Fermi interaction still doubt. 
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looked several ways. Its interpretation terms functional potential 
well has already been mentioned. Its similarity some the equations 
current field theory may also noted—including the nonlinearity; the 
other hand differs from these equations not being integral equation 
well. Finally, resembles especially the equations studied Reference 
(8); fact, the simplest nonlocal generalization those equations. 

The present investigation was intended primarily methodological 
interest. was showing that equation [1] could treated the 
same phase space methods were used earlier for the simpler equations 


described Reference (8). 
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